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We study the input-output response of a streamwise constant projection of the Navier-Stokes equations 
for plane Couette flow, the so-called 2D/3C model. Study of a streamwise constant model is motivated 
by numerical and experimental observations that suggest the prevalence and importance of streamwise and 
quasi-streamwisc elongated structures. Periodic spanwise/ wall- normal (z-y) plane stream functions are used 
as input to develop a forced 2D/3C streamwise velocity field that is qualitatively similar to a fully turbulent 
spatial field of DNS data. The input-output response associated with the 2D/3C nonlinear coupling is used 
to estimate the energy optimal spanwise wavelength over a range of Reynolds numbers. The results of the 
input-output analysis agree with previous studies of the linearized Navier-Stokes equations. The optimal 
energy corresponds to minimal nonlinear coupling. On the other hand, the nature of the forced 2D/3C 
streamwise velocity field provides evidence that the nonlinear coupling in the 2D/3C model is responsible 
for creating the well known characteristic "S" shaped turbulent velocity profile. This indicates that there is 
an important tradeoff between energy amplification, which is primarily linear and the seemingly nonlinear 
momentum transfer mechanism that produces a turbulent-like mean profile. 



I. INTRODUCTION 

We study the input-output response of the Navier- 
Stokes (NS) equations for plane Couette flow. Although 
this type of analysis can be performed for a variety of 
input/output combinations for the full nonlinear equa- 
tions, the mathematical complexity of such an endeavor 
makes it difficult to both obtain and interpret the results. 
Instead we perform the analysis in the simplified setting 
of a nonlinear streamwise constant projection of the NS 
equations. 

The choice of a streamwise constant model is moti- 
vated by studies of the linearized Navier-Stokes (LNS) 
equations, which show that streamwise constant features 
are the dominant mode shapes that develop under various 
perturbations about both the laminar— ~— and turbulent 
mean velocity^ profiles. In addition, streaks of stream- 
wise velocity naturally arise from the set of initial condi- 
tions that produce the largest energy growth 7 -^, namely 
streamwise vortices. Even in linearly unstable flows, 
studies have shown that the amplitude of streamwise 
constant structures can exceed that of the linearly un- 
stable modes^^. Bamieh and DahlcbJ, explicitly showed 
that streamwise constant perturbations produce energy 
growth on the order of i? 3 whereas disturbances with 

3 

streamwise variations produce growth on the order of R = . 

The prevalence of large-scale streamwise constant 
structures is also supported by direct numerical sim- 
ulation (DNS) and experiments. In Couette flow, 
DNS has long produced turbulent flows with very large 
scale streamwise and quasi-streamwise structures in the 
cor o 11 ' 12 . Experimental high Reynolds number studies 
have similarly identified large-scale streamwise coherence 
in other flow configurations^— . At high Reynolds num- 



bers (e.g., Re T > 7300), there is experimental evidence 
from pipes and turbulent boundary layers suggesting that 
these structures contain more energy than those nearer 
to the wal l 14 ' 16 ' 17 . The near-wall features, which also 
exhibit streamwise and quasi-streamwise alignment, are 
known to play a key role in energy production through 
the well studied "near-wall autonomous cycle"— ~—. 

The dominance of streamwise constant features was 
previously used to motivate the study of a streamwise 
constant model for plane Couette flow^. A stochastically 
forced version of this streamwise constant projection of 
NS was shown to reproduce important features of fully 
developed turbulence, including the shape of the turbu- 
lent velocity profile. Further, using Taylor's hypothesis, 
the same model also generated large scale streaky struc- 
tures, that closely resemble large-scale features in the 
core 23 Qj ven that maximum amplification of the LNS 
also occurs for the k x = modes (i.e., in a streamwise 
constant sense), the analysis of a streamwise constant 
model can be viewed as a study of the full system along 
the direction of maximum amplification. This approach 
allows us to understand the interaction between the well 
studied linear amplification mechanisms and additional 
effects due to the nonlinear coupling in the streamwise 
constant model. 

This paper is organized as follows. In the next sec- 
tion, wc describe the streamwise constant (so-called 
2D/3C) model and the idealized steady-state stream 
function model representing the cross-stream compo- 
nents of streamwise homogenous features. This stream 
function is used as input to a steady-state 2D / 3C stream- 
wise momentum equation. The solution corresponding to 
each stream-function input can be thought of as a forced 
solution of the respective streamwise deviation from the 
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FIG. 1. (a) Contour plot of the first y harmonic (q = f ) for the stream function model. This model represents the streamwise 
constant streaks and vortices commonly observed in DNS and experiments, (b) The stream function computed based on the 

a;-averaged spanwise DNS velocity field, which was integrated to obtain the stream function, i.e., ip Xave (y, Z) = These 
plots are reproduced from22. 



laminar flow. These streamwise velocity fields are com- 
pared to a spatial field of DNS data obtained from the 
Kawamura group21 in order to verify the ability of the 
model to capture the relevant features of turbulent flow. 
We compute the spanwise/ wall- normal (z-y) plane forc- 
ing required to produce each of the stream functions de- 
scribed above and study the input-output response from 
this forcing input to the streamwise velocity's deviation 
from laminar. The optimal spanwise wavelength com- 
puted in this manner is consistent with linear studies. 
However, the nonlinearity in the model gives additional 
insight to the relationship between amplification and the 
turbulent velocity profile. In fact, this work demonstrates 
that there is an important tradeoff between linear ampli- 
fication mechanisms and the nonlinearity required to de- 
velop an appropriately shaped turbulent velocity profile. 
The paper concludes with a summary of our results and 
directions for future work. 



II. MODELS 

A. The 2D/3C Model 

The 2D/3C model for plane Couette flow discussed 
herein is obtained by setting the streamwise (x-direction) 
velocity derivatives in the full NS equations to zero^. 
This can be thought of as a projection of the NS into 
the streamwise constant space. The velocity field is then 
decomposed into components u = [U + u' sw , V + v' sw , W + 
w' sw }; where (U, V, W) with U = U(y) = y, V = W = 
is the laminar flow and (u' sw , v' sw , w' sw ) are respectively 
the streamwise, wall-normal and spanwise time depen- 



dent deviations from laminar in the streamwise constant 
sense. One can explicitly show that for Couette flow this 
2D/3C formulation also results in a system with zero 
streamwise pressure gradient. 

A stream function ip[y,z,t), such that 

, dtp , dip 

Vsw = ~dl' Wsw = ~~dy~ 

ensures that the resulting model satisfies the appropriate 
2D continuity equation. This yields 



du' 



dt 
dAtp 
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+ gis-. There is no slip or penetration at 
the wall and periodic boundary conditions are assumed 
for the spanwise direction. 

The Reynolds number employed for all computations 
described herein is R = R w = ^f^, where the U w is the 
velocity of the top plate, h is the channel height and v 
is the kinematic viscosity of the fluid. All distances and 
velocities are respectively normalized by h and U w . 



B. The Stream Function Model 

As a first step, we focus on the effect of large-scale 
streamwise elongated features in the core of a fully turbu- 
lent flow. We limit our study to cross-stream (i.e. wall- 
normal/spanwise plane) inputs because energy amplifi- 
cation of perturbations (forcing) from the wall-normal 
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FIG. 2. Contour plots of (a) u' aWss , from ip as (y,z) = 0.00675 sin 2 (7ry) cos (y^2 z ) anc ^ C 3 ) "L t the streamwise velocity 
component of the rr-averaged DNS data. All plots correspond to R = 3000 and have the same contour levels. These plots are 
reproduced from 2 ^. 



and spanwise directions have been shown to scale as R 3 
whereas all of the other input-output combinations admit 
only R scaling^. We focus on the effect of these cross- 
stream inputs on the streamwise component of the flow. 

We are interested in developing a simple analytic 
model for the steady-state stream functions tpssil/t z ) that 
will define our inputs. This will lead to computational 
tractability and better lends itself to analytical studies. 
In Barkley and Tuckerman 2 ^ it was shown that laminar- 
turbulent flow patterns in plane Couette flow can be re- 
produced using a cross-stream stream function of the 
form ip(y, z) = -0o(y) + cos(fc z z) + ip 2 (y) sin(fc z z). 

We use this study as guidance but set the zeroth-order 
term to zero because a nonzero ipo produces a nonzero- 
mean spanwise how w' ss , which is not representative of 
the velocity field we are interested in studying. This leads 
to the following simple doubly harmonic function 



[ ^-z 



(2) 



.22,27 



which obeys the boundary condition: 

In order to estimate the values for X z , e and q we av- 
eraged a full spatial field of DNS data^ along the x- 
componcnt to obtain an approximation of a streamwise 
constant flow field. In the sequel, we refer to the this 
streamwise averaged DNS field as the x-averaged DNS 
data. A full discussion of this velocity field and its use as 
an approximation for streamwise constant data is given 
i n 22 i 27 . The x-averaged flow fields were used to ensure 
that our ip ss model as well as the corresponding wall- 
normal and spanwise velocities, have the correct features. 

Figure 1 (b) shows tp = — dW Qy Vi computed by integrat- 
ing the re-averaged spanwise velocity (w' Xave {y, z)) from 
the DNS data at R — 3000. The ip ss model correspond- 
ing to equation @ with q = 1, s = 0.00675 selected to 



match the magnitude of the integrated w' x and v' x 
DNS fields and A z w 1.8 selected to match the DNS' 
fundamental spanwise wavelength is provided in Figure 
1(a) Comparison of figures 1(b) and 1(a) indicate that 



the tpss model shows good agreement with the DNS data 
approximation in the region of highest signal. 



III. FORCED SOLUTIONS 

The steady-state solutions of ([Ta|) corresponding to 
steady-state stream functions of the form ((2|) are of in- 
terest for two reasons. First, they allow us to isolate 
the nonlinear streamwise velocity equation in order to 
demonstrate that its nonlinear coupling filters an appro- 
priately constructed tp ss (y,z) towards the expected "S" 
shape of the turbulent velocity profile. It also gives in- 
sight into mechanisms that create the momentum (en- 
ergy) transfer which generates this blunted profile. 

The forced solutions shown herein are presented solely 
to demonstrate that the 2D/3C model ([TJ is representa- 
tive of certain aspects of turbulent behavior and as such, 
an amplification study based on this model is of interest. 
Gayme et al^ 2 - and Gayme^i provide a detailed explo- 
ration of the extent to which the 2D/3C equations ([1]) 
can be used as a model for turbulent behavior. For all 
of the results presented in this section we solved for the 
forced solution u' sw (y, z) using both a least-squares ap- 
proach and itcratively using an explicit Eulcr method for 
comparison. The initial studies were carried out using the 
same grid resolution as in the DNS data described i»2A, 
(i.e. using a 96 x 512 grid on the y-z plane). We then 
reduced the resolution to a y-z plane grid of 48 x 100. We 
found negligible differences in the results between these 
two grid sizes. In the sequel, we only report the results 
for the 48 x 100 grid and the explicit Euler iterative so- 
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lution. 



A. The Velocity Field 



Figure 2(a) shows a contour plot of u' sw (y, z) result- 
ing from a stream function ip ss with the same param- 
eters as in Figure 1(a) (i.e. q = 1, s = 0.00675 and 
X z = 1.8, at R = 3000). This figure shows that our 
computed streamwise deviation from the laminar veloc- 
ity has features consistent with the difference between a 
laminar and turbulent velocity field. For example, the 
velocity gradients near the walls that are associated with 
an "S" shap ed (or blunted) turbulent velocity profile. 
Figure 2(a) also shows good qualitative agree ment with 
the s-averaged DNS data (shown in Figure |2(b)| with 
the same contour levels as 2(a)). This DNS data has 



previously been shown to have a spanwise mean veloc- 
ity corresponding to the full turbulent velocity profile^. 
The steady-state streamwise velocity deviation from lam- 
inar u' sw (y, z) has similar structural features to the DNS 
data in that both have near- wall minimum and maximum 
peaks that are out of spanwise phase with one another 
top-to-bottom. There is however, more variation in the 
deviation from laminar across the span as compared to 
the DNS field. 



B. Mean Profile 

The effective energy redistribution through the forced 
streamwise velocity deviations is further investigated 
through a spanwise mean over the streamwise velocity de- 
viation from laminar. This is carried out at five perturba- 
tion amplitudes (0.000675 < e < 0.02), all at R = 3000. 
Figure [3] shows averages across the span of u' sw (y,z) for 
these five e values along with a similar average of the 
x-averaged streamwise velocity field of the DNS. The 
spanwise average of the DNS, has been validated against 
other results in the literature by Tsukahara et al^. The 
use of ip ss from @ as an input to a the steady-state (fTa|) 
produces streamwise velocity profiles whose shapes are 
consistent with the a;-averaged DNS data. The peaks 
are, however located at different wall-normal positions. 

The fact that an exact match (with DNS data) for the 
wall-normal peak position is not obtained is not unex- 
pected given the simplicity of the wall-normal variation 
in the steady-state model, as well as the streamwise con- 
stant and steady-state assumptions. Clearly the full tur- 
bulent field is neither streamwise constant nor steady- 
state. The main point of presenting the velocity field 
arising from the forced steady-state model is to illustrate 
its effectiveness in reproducing the momentum redistri- 
bution associated with the change in the velocity profile 
from laminar to turbulent. This means that (jTJ) pro- 
vides more information about the turbulent velocity field 
than a linear model, (which cannot produce the change 
in mean velocity profile between laminar and turbulent 
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FIG. 3. Variation of the 2D/3C (streamwise constant) devi- 
ation from laminar, u' SWss with perturbation amplitude (e); 



based on input tp aa (y,z) 
is reproduced from22. 
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flows). Therefore studying input-output amplification in 
this model may provide us with some additional insight 
compared to the traditional analysis performed using lin- 
ear models. 

The simple steady-state model described herein rea- 
sonably predicts the essence of the mean behavior at the 
expense of losing some of the smaller scale details. For ex- 
ample, an exact characterization of the wall-normal vari- 
ation and, of course, the small-scale turbulent velocity 
fluctuations are not captured in this analysis. These re- 
sults suggest that the phenomenon that is responsible 
for blunting of the velocity profile in the mean sense is 
a direct consequence of the interaction between rolling 
motions caused by the y-z stream function and the lami- 
nar profile. In other words, this study provides strong 
evidence that the nonlincarity needed to generate the 
turbulent velocity profile is dominated by the nonlinear 
terms that are present in the u' sw (y, z,t) evolution equa- 
tion (TTal). 



IV. INPUT-OUTPUT AMPLIFICATION 

A. Energy Amplification 

In order to discuss input-output amplification, it is use- 
ful to determine the forcing required to produce a steady- 
state ipss- This is accomplished by solving a forced ver- 
sion of the steady-state ip evolution equation in (|lbj) for 
some forcing term r] S s{y,z). The linearized version of 
this forcing, which by abuse of notation we also refer to 
as r] ss (y, z), is given by 



Vss{y,z) = -—A 2 tp s 



(3) 



This rj S s can be viewed as the deterministic forcing re- 
quired to produce a particular ip ss ■ In the sequel, we use 
the linear rj ss of Equation ([3]) for all of the computations. 
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For a complete discussion of the use of a linear ijj equation 
see^. 

The input-output response can now be studied through 
an amplification factor of the form 



r M = 



(4) 



T ss is a nonlinear analog of the type L2~X>o-Li induced 
norm that has been used to study the optimal response of 
the LNS with harmonic input /forcing, see for example 2 ^. 
The energy in (j4]) is defined in terms of the squared 2- 
norm. For each 2-dimensional component p(y, z) this 
quantity is 



1011 



zi JO 



P(y, z ) 2 dy dz 



AzAy 



JV.-l Ny-l 



(5) 



J2 E Aw.**)- 



fc=l i=l 



where Ay = yi~ y\ and Az = Z2 — Zi respectively rep- 
resent the space between the z and y grid points. 



B. Reynolds Number Scaling 

The scaling of u' sw with R for a particular e is un- 
clear from equations (jlaj) and ([2]). An empirical rela- 
tionship was computed using the stream function model 
@ with q = 1 and e = 0.001 for four different values 
of i?: 3000, 6000, 10000 and 12000. Figure |4(a)| shows 
that \\u' sw || 2 appears to scale as a function of y/R at the 
higher wavenumbers (k z > 2), for the R values selected. 
The energy peaks also seem to collapse under this scaling 
for the higher Reynolds numbers we considered. 

The scaling of T ss can be estimated by combining this 
vR scaling of ||u' su , ss || 2 with the scaling of r) ss (y, z), as 

seen in ([3]). Thus, T ss should scale as a function of R^ . 
Figure |4(b)| which shows ^§f^ for the same R's and e 
indicates that the amplification factor data does collapse 
well under the i? 5//2 scaling, especially at the higher wave 
numbers. However, the R = 3000 data peak does not 
seem to follow this relationship. This discrepancy can be 
explained by looking at previous linear studies. 

The scaling of the input-output amplification for 
streamwise constant disturbances of the LNS equations 
was expressed as fi(k z )R + f2(k z )R 3 by Bamieh and 
Dahleh 2 -. Furthermore, they reported that the form of 
fi(k z ) and f2(k z ) means that the linear term dominate 
at low Reynolds numbers 2 -. The difference in scaling over 
different Reynolds number ranges was confirmed by a 
low Reynolds number linear study of Poiseuille flow that 
showed energy amplification at k x = scales with R? for 
the range 800 < R < 5000 and R 3 for larger Reynolds 
numbers 2 ^. In that study, R was normalized on half chan- 
nel height S. The equivalent normalization would make 



our Reynolds number range 750 < R c < 3000. Opti- 
mal amplification studies based on initial conditions also 
support R scaling at low Reynolds numbers 2 ^. Based on 
these studies, both the fact that our scaling is less than 
R 3 and that we found a lower T ss peak value for R = 3000 
(corresponding to R c = 750) are not unreasonable given 
the low Reynolds numbers we are employing. 



C. Optimal Spanwise Spacing 



Figure 4(a) indicates that \\u' sl 



' increases with k z 
until it reaches a maximum value and then levels off. 
We can similarly find a relationship between k z and T ss 
by substituting the expression for ^> ss from into the 
linearized noise equation ([3]) 



Vss{y, z ) = ^ cos(k z z)\ (8<zV 4 + 2q 2 k 2 z n 2 ) + 



4g 2 /c 2 7r 2 - (fc 2 + 4<2 2 7r 2 ) 2 sin 2 (qwy) 



(6) 



Equation ([6]) illustrates that the noise scales with k z and 
q 4 . So, the forcing energy ||?7 S s|| 2 monotonically increases 
with k z while ||m' su , ss || 2 peaks and then levels off. Thus, 
even though larger k z is associated with higher forcing 
the corresponding amplification factor does not continue 
to increase. There is an optimal k z that generates the 
most amplification: This is the dominant wavenumber 
corresponding to optimal spanwise spacing. In this sec- 
tion we explore how changes in e and R relate to the 
optimal spacing and discuss how our results compare to 
what has been previously reported in the literature. 

The peak values of T ss for the Reynolds numbers 
considered in Figure |4(b)| correspond to spanwise wave 
numbers of k z — 0.86, 1.0, 1.4 and 1.7; for R = 
12000, 10000, 6000 and 3000 respectively. This amounts 
to wavelengths of 7.3/i, 6.1/i, 4.6ft and 3.7ft. To determine 
if these are optimal values over the entire parameter set 
we need to determine the relationship between T ss and 
amplitude s. 

Figure 4(c) shows T ss for an amplitude range of 
0.000675 < e < 0.005 all with R = 3000. Over most of 
the range both r ss and the optimal spanwise wavenumber 
monotonically decrease with increasing e, although there 
appears to be a collapse at the minimal wavelengths. 
Therefore, the peak T ss over all the amplitudes we se- 
lected (i.e. the optimal k z ) occurs at the lowest e. Figure 
|4(d)| shows that continuing to reduce e results in conver- 
gence to an optimal wavenumber of k z = 2.06, which 
corresponds to X z = 3.05ft, for all e < 0.0001. We ob- 
tain the same optimal wavelength when we repeat this 
procedure for R = 6000, 10000 and 12000. 

Much of the literature, (e.g.,— iA£) related to optimal 
spanwise spacing has shown optimal spanwise wave num- 
bers k z £ [2.8,4]. However, in many cases these studies 
were aimed at determining spanwise spacing in the near 
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(c) (d) 

FIG. 4. (a) The streamwise energy scales as -/R. (b) The amplification factor r ss scales as . The optimal spanwise 
wavenumber occurs at the maximum F ss for each R. (c, d) T 3S for different values of e all at R — 3000. Both T 3S and the 
optimal spanwise wavenumber monotonically decrease with increasing e. 



wall (inner scaled) region. Recent Poiseuille flow stud- 
ies using the LNS linearized about a turbulent velocity 
profile, where an eddy viscosity is used to maintain the 
profile, found that at high Reynolds numbers there are 
actually two peaks in the optimal energy growth curves, 
one scaling in inner units and the other in outer units^^. 
The outer unit peak appears to correspond to the large- 
scale structures that have a spanwise spacing of approx- 
imately \ z G [2,5.2]<5. 

The only Couette flow study to look at both inner and 
outer unit scalings reported results at R = 300021. At 
this low Reynolds number there is little to no scale sep- 
aration between the peaks. They studied several input- 
output response functions and found that the optimal 
response to harmonic forcing occurs when X z = 3.85h. 
This harmonic forcing study is more closely related to our 
analysis than the initial condition-based studies reported 
in most of the other work. Couette flow DNS^ and ex- 
perimental studies^ have observed large-scale (outer re- 
gion) feature spacing in the range k z = [2,2.55]/i. Thus, 
our results (A 2 = 3.05/t.) lie right in between those of the 
analytical Couette flow study^i and the flow field obser- 
vations. A constant optimal wavelength across Reynolds 



numbers is also consistent with previous studies using a 
linear model with an eddy viscosity based turbulent ve- 
locity profile£i2£. 



D. Mean Velocity Profile versus Optimal k z 



Figure 5(a) shows the steady-state mean velocity pro- 
file computed using tp ss models ^ with amplitudes in 
the range 0.001 < e < 0.00675 at their corresponding 
optimal values of k z along with the DNS data. All plots 
correspond to R = 3000. The velocity gradients at the 
wall increase with e. The fit in the center of the channel 
also approaches the DNS data as e increases, although at 
e = 0.005 and above the curve overshoots the DNS. There 
is no amplitude that exactly matches the DNS data and 
the fit is especially bad in the near-wall region. As previ- 
ously discussed, this is because the assumptions inherent 
in the 2D/3C model neglect the smaller scale activity 
that dominates in the near-wall region. In the core, the 
mean velocity curves, ^° ss , corresponding to e = 0.005 
and e = 0.00675 respectively cross the DNS curve at a 
y + w 30 and y + ss 27, based on the DNS viscous units. 



(a) 



(b) 



FIG. 5. (a)The mean velocity profile of DNS data along with one computed from the steady-state (|la[) for a ip ss model ([2]) 
with 5 = 1, over a range of e with k z corresponding to the peak Tss for each e considered, (b) The mean velocity profile for 
e — 0.00675 at a number of different k z values compared with DNS data. The data in both (a) and (b) correspond to R = 3000. 




FIG. 6. The velocity gradient at the wall continues to increase 
while both T ss and the energy ||usii. ss || 2 peak and then drop 
off. The solid black line represents the peak F ss for each e. 



The maximum overshoot in the core (defined by y + > 30 
in DNS viscous units) is 3.6% and 6.2%, respectively for 
e = 0.005 and e = 0.00675. This is remarkably good for 
such a simplified steady-state model. 

The fact that the e = 0.005 and e = 0.00675 mean ve- 
locity profiles show the best agreement (most blunting) 
with the DNS data is not consistent with the fact that 
maximum amplification occurs at the smallest ampli- 
tudes (i.e. e < 0.0001). In order to study this further we 
looked at different k z values corresponding to a ip ss model 
amplitude of e = 0.005. Figure 5(b) shows the mean ve- 
locity profile of the DNS along with mean velocities for 
e = 0.005 at the maximum T ss (optimal wavenumbcr 
k z = = 0.69), at k z = 1.8 and at k z = 3.4. This last 
value coincides with X z — 1.8, i.e., the value correspond- 
ing to the dominant wavenumber of the x-averaged DNS 
data— and the results discussed in Section IIIII Again, 
while by definition the amplitude of T ss is larger for the 
optimal wavenumber k z = 0.6825, the velocity profile has 
larger velocity gradients at the wall for higher values of 
k z . This continued increase in shear stress at the wall as 
both k z and e increase is better seen in Figure [5] which 



depicts the energy in u' sa versus the shear stress at the 
wall. 

Small tp ss model ((2|) amplitudes, e's, correspond to 
lesser nonlinear coupling between the equations. On the 
other hand, at higher e the nonlinear terms have a larger 
magnitude because e directly multiplies each of the non- 
linear terms. As e decreases the energy amplification 
increases but the velocity profile blunting decreases (i.e., 
the profiles become increasingly laminar-like). This be- 
havior can be interpreted as follows; the amplification is 
dominated by linear mechanisms, whereas the blunting 
comes from nonlinear interactions. Moreover, there is 
some tradeoff between energy amplification and the cre- 
ation of a turbulent-like blunted mean velocity profile. 
The observation that blunting continues to increase with 
wavenumbers beyond the energy optimal wavenumber as 
depicted in Figure |6] appears to indicate that the exact 
relationship depends on the spanwisc wavenumbcr. This 
type of dependence is consistent with linear amplification 
theory. Further understanding of this tradeoff may pro- 
vide important insight into the mechanisms associated 
with both transition and fully turbulent flow. 



V. CONCLUSIONS AND FUTURE WORK 

A simple cross-stream model of large-scale stream- 
wise elongated structures nonlinearly coupled through a 
steady-state 2D/3C streamwise momentum equation al- 
lows us to isolate important mechanisms involved in de- 
termining the shape of the turbulent velocity profile. The 
momentum redistribution that produces features consis- 
tent with the mean characteristics of fully developed tur- 
bulence appear to be directly related to the 2D/3C non- 
linear coupling in the streamwise velocity evolution equa- 
tion. The steady-state 2D/3C model produces a blunted, 
turbulcnt-likc profile using very simple stream functions. 
This behavior appears to be robust to small changes in 
the stream function model. This repeatability suggests 
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a preference for redistribution of momentum along the 
wall-normal direction. Further understanding of the un- 
derlying dynamics of this mechanism may provide insight 
into the transition problem and allow better design of 
turbulent suppression flow control algorithms. 

An input-output analysis in this framework not only 
provides results consistent with previous studies but also 
illuminates an interesting interaction between energy am- 
plification and the increased velocity gradient at the wall 
associated with the turbulent profile. Essentially, al- 
though the input-output amplification monotonically de- 
creases with increasing forcing amplitude, the velocity 
profiles become increasingly more blunted. Thus, there is 
likely a tradeoff between the linear amplification mecha- 
nisms and nonlinear blunting mechanisms that determine 
the nature of the turbulence-like phenomena modeled by 
(fTj). This tradeoff may have important implications for 
flow control techniques that target skin friction or the 
mean profile. 

A natural extension of this work would be to refine our 
stream-function model through adding additional wall- 
normal and spanwise harmonics to the existing form of 
Equation ([2]). Developing an entirely new model using 
models of real sources of flow disturbances, as discussed 
ir*22, may also provide guidance in determining a better 



ip ss as a forced steady-state solution of (|Tb 



VI. ACKNOWLEDGEMENTS 

The authors would like to thank H. Kawamura and T. 
Tsukahara for providing us with their DNS data. This 
research is partially supported by AFOSR (FA9550-08-1- 
0043) and B.J.M. gratefully acknowledges support from 
NSF-CAREER award number 0747672 (program man- 
agers W. W. Schultz and H. H. Winter). 

1 B. F. Farrell and P. J. Ioannou, "Stochastic forcing of the lin- 
earized Navier-Stokes equations," Phys. of Fluids A, 5, 2600- 
2609 (1993). 

2 B. Barnich and M. Dahleh, "Energy amplification in channel 
flows with stochastic excitation," Phys. of Fluids, 13, 3258-3269 
(2001). 

3 M. R. Jovanovic and B. Bamieh, "The spatio-temporal impulse 
response of the linearized Navier-Stokes equations," in Proc. 
American Control Conf. (Arlington, VA, USA, 2001) pp. 1948- 
1953. 

4 M. R. Jovanovic and B. Bamieh, "Componentwise energy ampli- 
fication in channel flows," J. Fluid Mcch., 534, 145-183 (2005). 

5 J. C. del Alamo and J. Jimenez, "Linear energy amplification in 
turbulent channels," J. Fluid Mech., 559, 205-213 (2006). 

6 B. J. McKeon and A. S. Sharma, "A critical-layer framework for 
turbulent pipe flow," J. Fluid Mech., 658, 336-382 (2010). 

7 K. M. Butler and B. F. Farrell, "Three-dimensional optimal per- 
turbations in viscous shear flow," Phys. of Fluids A, 4, 1637-1650 
(1992). 

8 B. F. Farrell and P. J. Ioannou, "Optimal excitation of three- 
dimensiogal perturbations in viscous constant shear flow," Phys. 
of Fluids A, 5, 1390-1400 (1993). 

9 M. R. Jovanovic and B. Bamieh, "Unstable modes versus non- 
normal modes in supercritical channel flows," in Proc. American 
Control Conf. (Boston, MA, USA, 2004) pp. 2245-2250. 



10 L. H. Gustavsson, "Energy growth of three-dimensional distur- 
bances in plane Poiseuille flow," J. Fluid Mech., 224, 241-260 
(1991). 

11 M. J. Lee and J. Kim, "The structure of turbulence in a simulated 
plane Couette flow," Thin Solid Films, 1, 531-536 (1991). 

12 K. H. Bech, N. Tillmark, P. H. Alfredsson, and H. I. Ander- 
sson, "An investigation of turbulent plane Couette flow at low 
Reynolds numbers," J. Fluid Mech., 286, 291-325 (1995). 

13 K. J. Kim and R. J. Adrian, "Very large scale motion in the outer 
layer," Phys. of Fluids, 11, 417-422 (1999). 

14 J. F. Morrison, B. J. McKeon, W. Jiang, and A. J. Smits, "Scal- 
ing of the streamwise velocity component in turbulent pipe flow," 
J. Fluid Mech., 1508, 99-131 (2004). 

15 M. Guala, S. E. Hommema, and R. J. Adrian, "Large-scale and 
very-large-scale motions in turbulent pipe flow," J. Fluid Mech., 
554, 521-542 (2006). 

16 N. Hutchins and I. Marusic, "Evidence of very long meander- 
ing structures in the logarithmic region of turbulent boundary 
layers," J. Fluid Mech., 579, 1-28 (2007). 

17 N. Hutchins and I. Marusic, "Large-scale influences in near-wall 
turbulence," Phil. Trans. Royal Soc. London A, 365, 647-664 
(2007). 

18 F. Waleffe, "Proposal for a self-sustaining mechanism in shear 
flows," Tech. Rep. (Center for Turbulence Research, Stanford 
University/NASA Ames, 1990). 

19 J. M. Hamilton, J. Kim, and F. Waleffe, "Regeneration mecha- 
nisms of near-wall turbulence structures," J. Fluid Mech., 287, 
317-348 (1995). 

20 F. Waleffe, "On a self-sustaining process in shear flows," Phys. 
Fluids, 9, 883-900 (1997). 

21 J. Jimenez and A. Pinelli, "The autonomous cycle of near- wall 
turbulence," J. Fluid Mech., 389, 335-359 (1999). 

22 D. F. Gayme, B. J. McKeon, A. Papachristodoulou, B. Bamieh, 
and J. C. Doyle, "Streamwise constant model of turbulence in 
plane Couette flow," J. Fluid Mech., in press (2010). 

23 D. F. Gayme, B. McKeon, A. Papachristodoulou, and J. C. 
Doyle, "2D/3C model of large scale structures in turbulence in 
plane Couette flow," in Proc. of 6 th Int'l Symposium on Turbu- 
lence and Shear Flow Phenomenon (Seoul, Korea, 2009). 

24 T. Tsukahara, H. Kawamura, and K. Shingai, "DNS of turbulent 
Couette flow with emphasis on the large-scale structure in the 
core region," J. Turbulence, 7 (2006). 

25 K. M. Bobba, Robust flow stability: Theory, computations and 
experiments in near wall turbulence, Ph.D. thesis, California In- 
stitute of Technology, Pasadena, CA, USA (2004). 

26 D. Barkley and L. S. Tuckerman, "Mean flow of turbulent- 
laminar patterns in plane Couette flow," J. Fluid Mech., 576, 
109-137 (2007). 

27 D. F. Gayme, A Robust Control Approach to Understanding 
Nonlinear Mechanisms in Shear Flow Turbulence, Ph.D. thesis, 
California Institute of Technology, Pasadena, CA, USA (2010). 

28 P. J. Schmid and D. S. Henninson, Stability and Transition in 
Shear Flows (Applied Mathematical Sciences, Vol. 142, Springer- 
Verlag - New York, 2001). 

29 B. F. Farrell and P. J. Ioannou, "Variance maintained by stochas- 
tic forcing of non-normal dynamical systems associated with 
linearly stable shear flows," Phys. Rev. Letters, 72, 1188—1191 
(1994). 

30 G. Pujals, M. Garcia- Villalba, C. Cossu, and S. Dcpardon, "A 

note on optimal transient growth in turbulent channel flows," 

Phys. of Fluids, 21, 015109:1-015109:6 (2009). 
31 Y. Hwang and C. Cossu, "Amplification of coherent streaks in the 

turbulent Couette flow: an input-output analysis at low Reynolds 

number," J. Fluid Mech., 643, 333-348 (2010). 
32 J. Komminaho, A. Lundbladh, and A. V. Johansson, "Very large 

structures in plane turbulent Couette flow," J. Fluid Mech., 320, 

259-285 (1996). 

33 0. Kitoh and M. Umeki, "Experimental study on large-scale 
streak structure in the core region of turbulent plane Couette 
flow," Phys. of Fluids, 20, 025107-025111 (2008). 



